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SUMMARY 


The  Doppler-dif ference  optical  arrangement  is  nowadays  the  one  most  commonly 
used  for  fluid  flow  investigations  with  a  laser  anemometer.  In  this  review,  the 
essential  characteristics  of  the  Doppler-difference  signal  and  of  the  associated 
autocorrelation  function  are  considered  for  both  laminar  and  turbulent  flows. 
Experimental  conditions  under  which  it  is  known  that  the  autocorrelation  function 
can  be  analysed  to  yield  reliable  estimates  of  mean  velocity  and  turbulence 
intensity  are  specified. 


The  various  data-reduction  methods  which  have  been  proposed  are  classified 
and  briefly  reviewed.  It  is  shown  that  unknown  flows  having  relatively  high 
levels  of  turbulence  can  only  be  treated  successfully  if  there  are  very  many 
fringes  within  a  beam  diameter,  or  if  frequency-shifting  techniques  are  used  to 
achieve  an  equivalent  effect.  The  estimation  of  turbulence  intensity  in  low- 
turbulence  flows  also  presents  special  difficulties  and  a  procedure  which  uses 
information  already  available  in  the  transform  plane  is  described. 
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1  INTRODUCTION 

We  consider  in  this  review  the  problem  of  extracting  information  concerning 
the  mean  value  and  the  fluctuations  of  some  specified  component  of  velocity  of  a 
fluid  from  the  autocorrelation  function  of  the  intensity  of  the  radiation  scat¬ 
tered  from  a  localised  region,  the  measuring  volume,  by  small  particles  embedded 
in  the  fluid.  It  is  assumed  that  the  scatterer  follows  the  bulk  motion  of  the 
fluid  with  perfect  fidelity.  We  shall  be  concerned  with  the  output  of  an  ideal 
autocorrelator  operating  on  the  signals  obtained  with  a  Doppler-dif ference  opti¬ 
cal  system  which,  because  of  its  ease  of  alignment,  inherent  stability  and  the 
efficiency  of  incoherent  light-collection,  has  become  the  most  widely-used  laser 
velocimetry  arrangement  in  fluid-mechanical  applications.  Detailed  discussions 
of  the  autocorrelation  functions  encountered  in  other  forms  of  light-beating 
experiments,  such  as  those  using  homodyne  or  heterodyne  techniques,  can  be  found 
in  Ref  1 .  We  shall  not  consider  the  form  of  the  output  obtained  from  Doppler- 
difference  signals  with  a  clipping  correlator.  (The  formulae  given  here  should 
remain  valid,  however,  if  the  clipping  level  is  set  equal  to  zero  and  the  number 
of  counts  per  sample  time  is  very  much  less  than  one.) 

2  THE  SCATTERED  INTENSITY  AND  ITS  AUTOCORRELATION  FUNCTION  IN  LAMINAR  AND 

TURBULENT  FLOWS 

We  assume  that  the  light  source  is  a  linearly  polarised  laser  operating  in 
the  fundamental  (Gaussian)  TEMqq  mode.  Coherence  requirements  are  very  modest 
and  restriction  to  a  single  longitudinal  mode  of  operation  is  usually  unnecessary. 
We  assume  further  that  the  motion  of  the  scattering  particle  follows  that  of  the 
fluid  with  perfect  fidelity;  in  aerodynamic  applications  sub-micron  particles  are 
usually  necessary,  and  are  always  advisable  for  experiments  on  gaseous  flows. 

It  will  be  shown  later  that  a  small  angle  of  intersection  between  the  incident 
beams  is  an  important  feature  of  a  system  designed  for  use  on  a  completely 
general  turbulent  flow.  Under  these  conditions  it  is  a  reasonable  approximation 
to  assume  that  for  equal  incident  beam-powers,  the  scattered  light  intensities  in 
any  given  direction  are  equal.  For  larger  angles  of  intersection  the  depth  of 
modulation  of  the  Doppler-difference  signal  may  depend  significantly  on  the 
optical  geometry,  and  account  would  have  to  be  taken  of  the  effect  in  any  more 
general  theory. 

We  suppose  then  that  the  scatterer  radiates  with  equal  efficiency  from  the 
two  incident  beams.  We  allow  for  the  fact  that  these  may  not  be  perfectly 
balanced  by  introducing  the  quantity  p  ,  the  ratio  of  their  amplitudes.  We 
shall  assume,  however,  that  the  beams  have  the  same  characteristic  radius  r^  ; 
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this  is  the  radial  distance  at  which  the  amplitude  of  the  electric  field  has 

fallen  to  e  1  of  its  axial  value.  If  0  is  the  angle  of  intersection  of  the 

beams,  the  output  of  a  square-law  detector  (such  as  a  photomultiplier  tube)  on 

,  o 

which  the  scattered  radiation  falls  will  take  the  form 


I(x,y,z) 


where  (x,y,z)  are  the  coordinates  of  the  scatterer' s  position.  Iq  is  a  constant 

of  proportionality  which  includes  other  constants  such  as  the  collecting 

efficiency.  If  the  particle  is  moving  with  velocity  (v  ,v  ,v  )  and  at  time 

x  y  z 

t  =  0  is  at  some  general  point  ^q^q>Zq)  on  its  trajectory  (see  Fig  1),  then 


X  =  ?0  +  V 

y  =  n0  +  v 

z  =  ?0  +  v 


I(x,y,z)  is  now  a  function  of  time,  I(t)  .  The  frequency  of  the  cosine  term  in 
equation  (1)  is  seen  to  be  v  / s  ,  where  s  is  the  Doppler-dif ference  'fringe' 
spacing,  defined  in  the  usual  way: 

s  =  \j^2  sin  -|j 

I ( t )  represents  the  analogue  signal  which  would  be  obtained  if  the  scat¬ 
tered  radiation  weie  sufficiently  intense.  At  very  low  signal-levels  it  repre¬ 
sents  instead  the  probability  of  photon  detection  by  the  photocathode. 

The  output  of  an  autocorrelator  operating  on  such  a  signal  will  be  propor¬ 
tional  to  the  function  0 

00 

vO 

oo 

G(x)  =  J  I(t)I(t  +  T)dt  • 
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This  expression,  relating  to  a  laminar  flow,  has  been  evaluated  for  a  general 
optical  geometry  in  Ref  2.  The  result  is  complicated,  but  simplified  forms 
yielding  more  insight  into  its  physical  significance  have  been  obtained  for  par- 

3 

ticular  cases  and  are  presented  elsewhere  in  these  proceedings  . 

A  formula  applicable  to  three-dimensional  turbulence  was  also  derived  in 
Ref  2  by  making  the  simplifying  assumption  that  the  beam  intersection  angle  is 
very  small,  so  that  the  measuring  volume  becomes  approximately  cylindrical. 
Except  in  the  case  of  fully  forward  or  fully  backward  scattering,  the  detector 
can  usually  be  arranged  to  provide  a  well-defined  cut-off  in  the  direction  of 
the  beams.  We  summarise  here  the  main  feature  of  the  subsequent  argument. 

Equation  (1)  can  be  rewritten  as 


One  usually  possesses  some  a  priori  knowledge  of  the  mean  flow  direction 
and  we  will  assume  that  the  beams  are  oriented  in  such  a  direction  that  transits 
for  which  vx  is  large  are  relatively  uncommon.  The  angle  of  intersection  is 
now  supposed  so  small  that  for  the  great  majority  of  transits  the  quantity 
x  sin  0/2  is  negligible  compared  with  y  cos  0/2.  ^As  an  example,  suppose  the 
fringe  size  is  10  pm  and  A  =  0. 5  pm.  Then  0  ^  A/s  —  0.05  ^  3°  ,  and 
sin  0/2  =  0.025  .  If  the  detector  aperture  is  adjusted  to  make  x^^  ■  Ar^  , 
say,  we  have 


0.1  r„ 


Then  a  good  approximation  to  equation  (2)  is 


w 


1 
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f(t)  = 


exp 


2,2  2. 

—  (y  +  z  ) 

r0 


.  .  2  ^  -  2ir 

1  +  p  +  2p  cos  —  y 
s 


=  (1  +  p2) 

exp 

2  /  2  ^  2. 
(y  +  z  ) 

ro 

/,  .  2tt  \ 

11  +  m  cos  —  yl  , 


(3) 


where  m  is  the  Michelson  fringe  visibility: 


m 


2p 


1  +  p‘ 


Further  simplification  is  possible.  Suppose  the  point  of  closest  approach 
to  the  x-axis  has  coordinates  (x^.y^.z^)  .  The  vector  joining  this  point  to  the 
x-axis  is  normal  to  the  particle  trajectory  and  lies  in  the  direction 


Hence 


y„v  + 
J0  y 


z_v 
0  z 


0 


and  we  find  from  equation  (3) 


f(t)  =  (1  +  P2) 

exp 

2  t  2  .  2. 

%  +  z0> 
r 

0 

exp 


-  2 


2v 

1  +  m  cos  —  (yn  +  v  t) 
s  w0  y 


(4) 


The  autocorrelation  function  of  this  expression  has  been  shown  in  Ref  2  to 
be,  to  a  very  close  approximation  , 


G(t)  =  -y  rQ( 1  +  p2)2 


-Pi -4  (y0  +  Z0} 

ro 


(v  +  v  ) 
y  z 


2x-4 


x  exp 


,  2^2. 
vy  *  V  2 
l~2  t 

ro 


1  +  J  in  cos 


(tv) 


(5) 


f(t)  and  G(x)  represent  the  function  of  time  and  the  autocorrelation  function 
which  would  be  obtained  under  conditions  of  laminar  flow,  or  from  a  single 
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strong  scatterer.  We  remark  here  that  they  contain  similar  exponential  multiply¬ 
ing  factors,  functions  of  t  and  x  respectively,  so  that  high-pass  filtering 
of  either  f(t)  or  G(x)  in  order  to  obtain  an  estimate  of  v^  from  the  cosine 
term  would  introduce  similar  errors. 


The  notation  will  be  changed  at  this  point  to  conform  with  a  common  aero¬ 
dynamic  convention.  The  velocity  component  normal  to  the  fringe  planes,  v^  , 
will  be  denoted  by  u  and  the  component  perpendicular  to  the  plane  of  the 


beams,  v  ,  by 


The  third  orthogonal  component  will  be  denoted  by  w  . 


Then  equation  (5)  becomes 


G(x) 


+  v^ 


exp 


,  2  A  2.  2 

(u  +  V  ) X 


i  2  2irux 
4m  cos  - 

*■  a 


)  ■ 


(6) 


{ 


i 


where  the  various  constants  have  been  included  under  the  symbol  a^  .  If  G(x) 

represents  a  summation  over  a  number  of  transits,  a.  will  be  understood  to 

2  2  ° 

involve  an  average  value  of  the  quantity  (yQ  +  z^) ,  in  an  obvious  sense. 

In  turbulent  flow,  successive  particles  cross  the  measuring  volume  gener¬ 
ally  with  different  velocities,  and  the  resulting  autocorrelation  function  will 
be  an  integral  of  the  laminar  form  for  G(x)  of  (6)  over  the  probability 

density  function  p(u,v,w)  which  characterises  the  turbulence.  The  fact  that 

.  .  .  .  .  2 
the  scattering  volume  is  approximately  cylindrical  can  be  exploited  to  show 

that  the  composite  autocorrelation  function  -  the  output  of  the  correlator  in 

fully  turbulent  conditions  -  becomes 


00  OO  00 


H(x)  =  a 


1  /  /  /p(u’v’w) 


—00  “  00  “00 


exp 


,  2  2>  2 
(u  +  V  )x 


fl  +  $m^  cos  dudvdw 

{  s  / 


(7) 


where  a j  now  includes  constants  depending  on  the  numbers  of  particles,  the 
optical  geometry  and  the  duration  of  the  experiment.  In  deriving  (7)  it  was 
assumed  that  the  particles  are  completely  randomly  distributed  in  an  incompress¬ 
ible  fluid,  and  that  the  dimensions  of  the  measuring  volume  are  small  compared 
with  the  length  scale  of  the  turbulence.  The  dependence  of  the  instantaneous 

rate  at  which  transits  occur  on  the  fluid  velocity  leads  to  the  disappearance  in 

/2  2 

(u  +  v  )  appearing  in  (6). 
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ANALYSIS  OF  THE  AUTOCORRELATION  FUNCTION 


3. 1  Laminar  flows 


Here  equation  (6)  applies  and  data-reduction  can  proceed  by  straight¬ 
forward  curve- fitting.  Alternatively,  the  problem  can  be  regarded  as  a  special 
case  of  the  more  general  one  of  estimating  mean  velocity  in  turbulent  flow  and  a 
single  computer  programme  can  be  designed  to  cover  both  conditions.  Purely 
laminar  flows  are  in  any  case  rarely,  if  ever,  encountered  in  practice  and  an 
estimate  of  the  turbulence  level  present  is  usually  wanted. 

3.2  Turbulent  flows 

We  assume  that  fluctuations  may  be  occurring  in  all  three  components  and 
note  that  integration  over  the  w-component  in  (7)  can  be  carried  out  immed¬ 
iately,  since  the  kernel  of  the  integral  is  independent  of  w  . 

By  definition,  the  joint  probability  for  the  u-  and  v-components  is 


puv(u'v)  = 


00 

I  P(u,v,\ 


and  hence 


oo  CO 

f  f 

/  1  p  ,/u»v> 

exp 

J  J  uv 

-CO  —CO 

,  2  .  2v  2 

(u  +  V 


.  2  ?ttut\  ,  j 
+  im  cos  — - — I  dudv 


In  practice  there  will  always  be  a  nearly-constant  (dc)  level  added  to  the 
signal,  which  arises  mainly  from  background  radiation  present  during  the  experi¬ 
ment.  This  is  the  value  to  which  G(x)  will  tend  as  x  tends  to  infinity  and 
is  normally  included  in  the  data-reduction  procedures  as  another  unknown  variable 
whose  magnitude  has  to  be  determined. 

Equation  (8)  presents  us  with  an  impasse,  for  we  cannot  recover  the  two- 
dimensional  distribution  pyv  from  an  experimental  function  of  a  single  var¬ 
iable.  The  problem  centres  essentially  on  the  exponential  term,  since  this 
factor  couples  both  u-  and  v-components  into  the  measured  correlation  func¬ 
tion,  and  it  is  obvious  that  some  extra  information  is  required  if  progress  is 
to  be  made  at  all.  If  a  complete  mathematical  model  for  p^v  is  available, 
data-reduction  can  be  carried  out  by  curve-fitting  to  H(x)  .  This  will  not  be 
the  case,  however,  in  most  problems  of  any  real  intrinsic  interest  ( is  for  the 
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majority  of  experiments),  and  we  now  consider  special  cases,  involving  weaker 
assumptions,  for  which  useful  information  about  the  velocity  field  can  be 
recovered  from  the  experimental  function  H(x)  . 

3.2.1  puv  separable 
We  write 

p  (u,v)  =  p  (u)p  (v) 

ruv  *u  rv 

Note  that  this  implies  that  the  stress  term  u'v'  must  be  zero;  primes  denote 
the  zero-mean  fluctuating  parts  of  the  velocity  components  and  capital  letters 
their  mean  values: 


u  =  U  +  u’ 


v  =  V  +  v' 


For,  by  definition 


uv 


~  J  f  UVPUV ( u , v) dudv 


— OO  —00 


00  00 
=  J  upu(u)du  J 


(u)du  I  up  (v)dv 


But 


Hence 


=  UV 


uv  *  (U  +  u’ ) (V  +  v? ) 


UV  +  u’v' 


Equation  (8)  can  be  written 


H(t) 


OO  00  p 

ii  j  pv<v)  exp^_  dv  J  pu(u)  exp^_iiT"^  ^  + 


2hut  \  , 
cos  - Jdu 

8  J 


(9) 


Both  integrals  are  unknown  functions  of  t  ,  so  that  progress  cannot  be 
made  in  determining  the  characteristics  of  either  p^  or  pv  without  further 
assumptions. 
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(a)  u—  and  v-components  both  normally  distributed 

We  put 

(u  -  U)2 

2a2 
u 


pu(u) 


exp 


where  a  is  the  standard  deviation  of  the  fluctuations  in  u  ,  with  a  similar 
u 

expression  for  py(v)  . 

In  this  case  the  integrals  can  be  evaluated  explicitly,  and  we  find,  to  a 
very  good  approximation. 


H(t) 


exp 


(U2  ♦  V2)t2 


1  +  Jin  exp 


/  ,  2  2  2\ 


cos 


2itut 


for  values  of  t  such  that 


T2  «  ll° 

2  2  » 
< 


2 

ilo 

2  2 
a 

v 


(10) 


(ID 


(With  these  restrictions,  the  result  can  also  be  shown  to  hold  if  p  is  a 

UV  4 

bivariate  normal  distribution,  so  that  a  shear  stress  term  can  be  accommodated  . 
In  this  case,  puv  is  not  separable.) 

It  is  often  more  convenient  to  work  in  terms  of  the  turbulence  intensity 
n  ,  the  ratio  of  the  standard  deviation  to  the  mean,  rather  than  in  terms  of  the 
standard  deviation  itself.  For  example,  suppose  we  are  using  a  48-channel 
correlator  operating  at  50  ns  per  channel  in  an  experiment  on  a  supersonic  flow 
where  the  mean  velocity  is  400  ms'.  Then  Tmax  is  2.4  ys  and  for  a  beam  dia¬ 
meter  of  1  mm,  the  inequalities  would  require  that  both  nu  and  nv  are  less 
than  about  0.15  (15%)  or  so. 

(b)  Hermite  polynomial  expansion  for  p^(u):  v-component  normally 
distributed 

This  approximation"*  is  intended  for  applications  where  the  distribution  of 
one  velocity  component  only  may  be  highly  non-Gaussian  and  rests  on  the  assump¬ 
tion  that  Py(v)  is  a  Gaussian  and  that  p^(u)  can  be  written  in  the  form  of  a 
Gram-Charlier  expansion: 


o 

00 

VO 
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1  1 


P  (u,v) 
*uv 


2ir a  a 
u  v 


exp 


(u  ~  U)' 
2a2 


(v  -  V)‘ 


2a 


2L-1  ,  . 

I 

k=0  '  ' 


where  A^ 


and  He^  is  the  Hermite  polynomial  defined  by  the  relation 
He^(z)  =  (-)k  exp(£z2)  £exp(-  ^z2)J 

The  choice  of  L  in  the  upper  limit  (2L  -  1)  of  the  sum  will  depend  on 
the  amount  of  information  required  from  the  calculation.  The  first  few  poly¬ 
nomials  have  the  following  forms: 

HeQ(z)  =  1  ,  He | (z)  =  2z  ,  He2(z)  =  4z2  -  2  , 

He^Cz)  =  8z2  -  12z  ,  He^(z)  =  16z^  -  48z2  +  12 

The  double  integral  can  again  be  evaluated  explicitly  and,  using  again  the 
approximations  of  (11),  takes  the  form 


H(x) 


(12) 


The  A^  are  now  determined  by  curve-fitting.  Fig  2  of  Ref  5  illustrates  the 
way  in  which  a  much  better  approximation  can  be  obtained  to  a  bimodal  form  of 
p^(u)  using  the  Gram-Char lier  method  than  with  the  purely  Gaussian  model  of  the 
last  section. 


(c)  Eigenfunction  expansion  for  p^(u):  v-component  assumed  negligibly 

small 

Here  the  correlation  function  reduces  approximately  to  the  form 


12 


H(t)  -a j  I  pu(u) 

—co 

It  is  convenient  at  this  point  to  introduce  an  important  practical  consid¬ 
eration.  Because  the  kernel  of  the  integral  is  an  even  function  of  u  ,  any 
negative  values  of  u  which  may  occur  will  contribute  to  H(t)  in  exactly  the 
same  way  as  equal  positive  ones.  Hence  if  an  undistorted  form  of  p^(u)  is  to 
be  recovered,  the  u-component  must  always  be  of  the  same  sign;  ie  the  flow 
must  be  non-reversing.  (In  addition,  the  analysis  summarised  in  this  section 
rests  on  Hilbert-Schmidt  theory,  which  requires  that  the  kernel  be  symmetric  in 
u  and  t  .)  Let  us  suppose  then  that  u  is  always  positive,  so  that  the  lower 
limit  of  the  integral  can  be  set  equal  to  zero: 

00 

H(t)  =  a t  j  Pu(u) 

0 


l  (  i 

exp^-  — 2~j  ^1  +  im 


2  2ttut\ 


cos 


s 


) 


idu 


(13) 


The  procedure  now  is  to  try  to  find  a  set  of  functions  and  factors 

X^,  ,  eigenfunctions  and  eigenvalues,  which  satisfy  the  associated  Fredholm 
equation 


where  K(ut) 


00 

J  <(k  (u)K(ur)du  =  X^<|k(t) 
U 


Suppose  also  that  the 
arbitrary  function,  such  as 


form  a  complete  orthogonal  set.  Then  any 
p^(u)  ,  can  be  expanded  in  a  series  of  the  form 


Pu(u) 


OP 

z 


cj4>  j  (u) 


where  the  c.  are  unknown  coefficients  which  have  to  be  determined.  Substitut- 
J 

ing  this  expression  into  the  equation  for  H(x)  and  proceeding  formally,  we 
have 


13 


H(t)  =  a. 


00  00 

/  Zvi 

0  l_j=0 


(u) 


K(ut)du 


OO 

Zc.  A  ,d> . 

3  3  3 
j=0 


(t) 


Using  the  orthogonality  of  the  ^  ,  we  have 


OO  00  r  00 

J  H(x)4>k(T)dT  =  aj  J  I  CjX (t) 

0  0  j=0 


<j>k(T)dT 


=  a,  c,  A 


iVk 


so  that 


=  7X  !  H(THk(T)dT  ’ 

K  0 


and  finally 


3\ 

00 
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oo  a 

"u<u>  -  77  Z  rr  Vu>  / 

j=0  J  0 


H(x) <t> ^  (t)dT 
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A  complete  set  of  eigenfunctions  and  their  associated  eigenvalues  for  the 
Doppler-dif ference  kernel  lias  been  found.  They  are  discussed  in  detail  and 
their  application  to  experimental  data  described  in  another  contribution  to 
these  proceedings^. 

3.2.2  puv  not  separable 
Consider  again  equation  (8): 


H(t) 


j  p  (u,v) 

J  ruv 


exp 


,  2  A  2.  2 

(u  +  V  )t 


1  +  jm  cos 


2 _ 2ffutJ 


dudv 
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We  discuss  two  important  cases,  distinguished  by  the  behaviour  of  the 
exponential  term. 

(a)  Exponential  term  approximately  unity 

If  there  are  sufficient  fringes  across  the  beam  (rQ  >  s)  the  exponential 
term  will  be  almost  unity  over  the  observed  length  of  the  autocorrelation 
function  and  we  obtain  immediately  the  approximation 


H(t)  a„ 


/  Pu 


,  N  2ttut  , 

(u)  cos  -  du  +  a, 

s  l 


since 


Pu(u) 


00 

*J 


p  (u,v)dv 
ruv  ’ 


Note  that  it  is  unnecessary  to  make  assumptions  now  about  the  shape  of  the 
function  p^v  and  that  high  turbulence  levels  in  the  v-component  will  not 
affect  the  result.  Again  we  remark  that  the  correct  form  for  the  distribution 
Pu(u)  can  only  be  obtained  if  the  u-component  is  of  constant  sign.  Here  we 
assume  that,  if  necessary,  a  constant  known  bias  uffi  is  imposed  in  order  to 
make  u  +  uffi  positive  for  each  particle  transit.  This  can  be  achieved  by  intro¬ 
ducing  a  frequency-shifting  device  into  one  or  both  beams;  see  for  example  Ref  7 
for  a  review  of  available  methods.  We  can  then  without  loss  of  generality  again 
replace  the  lower  limit  of  the  integral  by  zero: 


1 


i 

i 


j 


H(x)  =  a.  /  p  (u)  cos  -  du  +  a  .  (15) 

l  J  u  s  1 

0 

Hence  Pu(u)  can  be  recovered  by  performing  a  simple  Fourier  cosine  transform 
on  the  data,  after  removal  of  the  additive  background  constant  .  The  latter 
can  easily  be  carried  out  automatically  by  the  data-reduction  programne  -  high 
accuracy  is  not  necessary. 

The  requirement  that  r^  >  s  implies  a  lower  limit  to  the  spatial  resolu¬ 
tion,  characterised  by  the  beam  diameter  2r^  ,  which  can  be  achieved  in  high¬ 
speed  flows.  With  a  minimum  sample  time  of  50  ns  and  beams  of  not  more  than 
I  mm  diameter,  the  condition  cannot  be  satisfied,  for  all  practical  purposes,  in 
most  supersonic  applications.  (In  nearly  laminar  conditions  frequency-shifting 
could  conceivably  be  used  to  bring  the 
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low-speed  flows  it  is  very  often  possible  to  arrange  the  experiment  appropriat¬ 
ely.  Note  also  that  if  the  speed  is  sufficiently  low,  frequency-shifting 

g 

methods  can  now  be  used  to  accommodate  very  high  levels  of  turbulence  ,  for 
which  flow  reversals  will  occur  from  time  to  time,  or  to  achieve  fine  spatial 
resolution. 


(b)  Exponential  term  varying  significantly:  narrow  distributions  of 

the  velocity  components 


These  conditions  are  met  with,  for  example,  in  supersonic  boundary  layers 
2 

at  low  Mach  numbers  .  In  this  case,  since  u  and  v  are  assumed  to  lie  within 
a  narrow  range,  the  exponential  term  in  (8)  varies  much  more  slowly  within  the 
integral,  for  a  given  value  of  r  ,  than  the  cosine  term  and  can  be  replaced  to 
a  good  approximation  by  the  value  it  takes  when  u  and  v  are  given  their  mean 
values,  U  and  V  .  It  can  then  be  taken  outside  the  integral  and  the  integra¬ 
tion  over  the  v-component  again  carried  out,  yielding  the  approximate 
expression 


H(x) 


- 

exp 

(U2  +  V2)x2 

2 

ro 

l  pu(u)^l  +  Jm2  cos  — du 


Since  the  turbulence  level  is  low  the  lower  limit  can  be  replaced  by  zero,  and 
finally  we  have 


H(x)  =  a 


exp 


(U2  +  V2)t2 


00 

l  PU(U)(]  + 


,  2 

£m  cos  — - — Jdu 


r> 


(16) 


If  the  mean  flow  direction  is  not  known,  two  fringe  systems  can  be  used 
to  obtain  estimates  of  U  and  V  .  For  narrow  nearly  symmetric  velocity  distri¬ 
butions,  the  positions  of  the  peaks  in  the  transform  plane  are  barely  affected 
by  the  exponential  term  and  a  linear  approximation  to  it  will  provide  good 
initial  values  for  U  and  V  .  The  computer  programme  can  be  designed  to 


determine  r^  from  these  estimates,  together  with  any  residual  background 


contribution  to  the  correlation  function.  The  magnitude  of  the  exponential  term 
can  now  be  calculated  for  each  value  of  t  and  a  modified  form  of  H(t) 
derived  from  which  the  damping  effect  of  the  exponential  has  been  eliminated. 

A  second  Fourier  transform  will  now  give  an  improved  estimate  for  p  (u)  from 


which  nu  can  be  determined  with  good  accuracy.  A  similar  calculation  can  be 


carried  out  to  obtain  py(v)  ,  and  hence  ny  ,  from  the  other  set  of  data. 
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3.3  Data-reduction  using  Fourier  transforms 

Clearly  the  most  satisfactory  experimental  technique  is  to  arrange  where- 
ever  possible  that  there  are  very  many  fringes  contained  within  the  beam  dia¬ 
meter,  since  analysis  can  then  proceed  with  the  minimum  of  a  priori  assumptions 
about  the  flow.  In  addition,  we  have  seen  in  section  3.2.2(a)  that  the  problem 
of  data-reduction  collapses  to  that  of  carrying  out  a  Fourier  transform  on  the 
experimental  data,  which  is  a  simple  and  rapid  computational  procedure.  We 
remark  here  that  the  fast  Fourier  transform  (FFT)  is  not  usually  the  appropriate 
technique,  since  it  lacks  flexibility  and,  at  low  turbulence  levels,  the  neces¬ 
sary  resolving  power.  However,  it  has  been  successfully  used  in  difficult 

9 

engineering  applications  ,  and  if  sufficient  signal  is  available  from  each 
scatterer,  so  that  the  correlator  can  be  operated  in  a  burst  mode,  the  FFT  is 
the  obligatory  method  for  processing  each  autocorrelation  function. 

We  have  also  seen  that  the  Fourier  transform  method  can  be  adapted  to 
cover  experiments  on  high-speed  flows  where  the  turbulence  levels  are  fairly 
low,  which  includes  many  of  the  types  of  flow  of  particular  interest  at  the 
Royal  Aircraft  Establishment.  For  these  reasons,  our  efforts  in  the  field  of 
data-reduction  have  been  concentrated  mainly  on  the  development  of  Fortran 
programmes  based  on  this  technique.  A  detailed  description  of  the  structure  of 
such  a  programme,  and  of  the  hardware  used  for  on-line  applications,  is  given 
in  another  paper  presented  at  this  conference^.  We  now  discuss  some  import¬ 
ant  characteristics  of  the  Fourier  transform  in  laminar  and  low- turbulence 
flows. 

3.3.1  Laminar  flows 

The  limited  number  of  channels  available  in  a  clipping  correlator  signifi¬ 
cantly  restricts  the  number  of  cycles  available  in  the  record  obtained  from  a 
laminar  flow.  The  rectangular  data-window  produces  an  effect  in  the  transform 
plane  which  is  exactly  equivalent  to  instrumental  broadening  in  classical 
spectroscopy.  The  delta-function  which  would  ideally  be  present  is  replaced  by 
a  function  having  a  sin  x/x  behaviour,  although  still  centred  at  the  same 
point  on  the  u-axis  .  Thus  an  accurate  estimation  of  the  velocity  can  be  made 
simply  by  finding  the  position  of  the  peak  of  the  Fourier  transform.  (We  ignore 
the  extremely  small  perturbation  which  will  arise  from  the  associated  sin  x/x 
function  centred  at  the  symnetric  point  on  the  negative  u-axis.)  We  shall  now 
show  that  the  sin  x/x  structure  can  also  be  used  to  decide  whether,  within  the 
limits  of  experimental  accuracy,  the  flow  is  laminar. 
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We  assume  that,  if  present,  the  beam  profile  term  (the  exponential  factor) 
has  been  removed  in  the  manner  described  in  section  3.2.2(b)  together  with  any 
residual  constant  term.  Hence  we  have  an  equation  of  the  form 

00 

G(t)  =  c  J  Pu(u)  cos  ^TTg—  du  (17) 

0 

where  c  is  some  constant.  The  inverse  transform  is,  by  the  Wiener-Khinchine 
theorem, 

00 

PU(U)  =  ~TTC  j  G(T)  C°S  "^T1  dT  ’  (,8) 

0 

Now  for  laminar  flow  at  velocity  Uq  , 

2TTU-T 

G(t)  =  cos  - 

s 


The  data  record  terminates  at  some  value  ,  say,  of  t  ,  so  that  the  com¬ 
puted  function  will  be,  writing  w  and  u>q  for  2ttu/s  and  2ttuq/s  respectively. 


p  (u)  =  — 

rU  7TC 


1  sinj(u.  -  o,0)Tm}  sinj(<o  +  «0)Tmf 


U  +  U)r 


At  u  =  uft  ,  the  ratio  of  the  first  term  to  the  second  is  at  least  2wr.r  . 
u  u  m 

But  the  time  per  cycle  is  2tt/wq  ,  so  that  the  number  of  cycles  contained  in  the 
record,  say  q  ,  is  and  the  ratio  is  at  least  4:rq  .  A  typical 

experimental  value  for  q  ,  using  a  48-channel  correlator,  would  be  about  8, 
so  that  the  first  term  would  be  more  than  100  times  the  maximum  value  of  the 
second  at  resonance.  Hence  around  u  *  Uq  we  take  as  a  good  approximation 


in  terms  of  q  . 


18 


The  first  zeros  on  either  side  of  the  peak  at  u  =  Uq  occur  when 

2irq(—  -  A  =  ±  ir 

VO  / 

ie  when 


Suppose  we  regard  the  curve  between  these  two  points  as  defining  a  velo¬ 
city  distribution.  Then  the  'turbulence  intensity'  (the  ratio  of  the  standard 
deviation  to  the  mean)  associated  with  such  a  distribution,  say  ip  ,  would  be 
given  by  the  formula 


a+ 

* 

JL 

J  (u  -  u0)2Pu(u)du 

a_ 

uo 

a+ 

j  Pu(u)du 

J 

where  ot+  ■  Uq| 

9 

a 

i 

ii 

c 

o 

Putting 

X 

■  -  o 

we  find 

Ip 

*  ^  ^  j 

(20) 


From  tables. 


4. 82q 


Hence  ip  will  depend  purely  on  the  number  of  cycles  in  the  record.  The 
progranme  can  easily  be  arranged  to  determine  this  quantity  as  a  by-product  of 
other  calculations,  and  it  can  be  compared  with  the  estimate  of  turbulence 
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intensity  generated  by  the  main  programme.  (This  searches  outwards  from  the 
peak  in  the  Fourier  transform  plane  and  also  uses  as  limits  for  the  integration 
the  points  at  which  the  distribution  first  crosses  the  u-axis  .)  If  the  two 
values  agree  closely,  it  can  be  concluded  that  the  flow  is  essentially  laminar; 
any  turbulence  present  in  the  flow  would  result  in  an  estimate  for  turbulence 
intensity  greater  than  iji  . 

3.3.2  Flows  with  low  turbulence  levels 

At  sufficiently  high  turbulence  levels,  the  autocorrelation  function  will 
have  decayed  away  effectively  to  zero  by  the  end  of  the  record;  for  a  48-channel 
correlator  operating  with  typical  settings  on  a  flow  having  a  fairly  symmetric 
unimodal  velocity  distribution,  the  level  below  which  truncation  effects  begin 
to  appear  is  about  4  or  5%.  Mean  velocity  estimates  are  very  little  affected 
at  such  low  levels,  but  the  accurate  determination  of  the  turbulence  intensity 
begins  to  pose  problems.  If  a  good  model  of  the  flow  is  available  for  curve¬ 
fitting  purposes,  it  should  be  possible  to  obtain  accurate  estimates  down  to  the 
laminar  condition;  for  example,  a  Gaussian  model  would  often  be  acceptable. 

Methods  for  extrapolation  of  the  autocorrelation  function,  using  analytic  con- 
11  12 

tinuation  or  maximum  entropy  techniques,  have  also  been  proposed.  However, 
they  would  all  require  a  separate  and  distinct  programme.  A  procedure  has  there¬ 
fore  been  devised,  based  on  an  extension  of  the  technique  described  in  the  last 
section,  which  uses  information  already  available  in  the  Fourier  transform  plane. 

The  basic  assumption  is  that  at  these  low  turbulence  levels  the  distribu¬ 
tions  are  effectively  Gaussian.  For  each  value  of  turbulence  intensity  f| 
over  the  desired  region  (typically  0  -  5%),  an  autocorrelation  function  having  a 
specified  length  of  q  cycles  can  be  computed;  this  represents  the  data  that 
would  be  obtained  in  an  ideal  experiment  on  such  a  flow.  Since  the  autocorre¬ 
lation  data  are  truncated,  the  Fourier  transform  will  consist  of  a  broadened 
central  lobe  with  oscillatory  wings.  From  this  central  lobe,  defined  by  the 
points  at  which  the  curve  first  crosses  the  u-axis  on  either  side  of  the  peak, 
an  apparent  turbulence  intensity  n  can  be  calculated.  It  is  found  as  a  matter 

of  experiment  that  if,  for  various  values  of  q  ,  the  quantity  qn  is  plotted 

3 

against  qnt  ,  the  curves  virtually  coincide,  and  a  single  look-up  table  can 
therefore  be  constructed  to  cover  a  range  of  values  of  the  product  qna  . 
Estimates  of  q  and  n  are  available  from  the  main  programme  and  the  refined 

3 

estimate  of  turbulence  intensity,  n  ,  can  be  made  part  of  the  output  data-set. 
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3.4  The  autocorrelation  function  for  a  single  transit  in  the  presence  of  noise 

It  has  been  remarked  previously  that  an  experimental  autocorrelation  func¬ 
tion  will  always  contain  a  reasonably  constant  additive  term,  usually  arising 
mainly  from  the  stray  radiation  present  in  the  experimental  environment.  If  it 
can  be  assumed  that  only  one  scatterer  has  contributed  to  the  data,  interesting 
information  about  the  circumstances  of  the  experiment  can  be  extracted  in  the 
following  way. 


Equation  (4)  can  be  written,  ignoring  the  phase  constant  (2u/ s)y ^  , 


Kt)  =  c. 


exp 


2(u2  +  v2)t2 


i  ^  2rrut\  A 

1  +  m  cos  - 1  +  c 


(21) 


where  Cj  =  1^(1  +  p  )  exp 


2(2  2\ 


and  c^  is  the  intensity  of  the  background  illumination,  which  is  assumed  to  be 
constant.  Then  if  the  experiment  lasts  for  a  time  T  ,  long  compared  with  the 
beam  transit  time  rQ//(.u2  +  v2)  ,  we  have 
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The  second  exponential  term  will  be  quite  negligible  in  any  real  experi¬ 
ment  and  we  can  write 


G(t) 


2  n  ro 


'i  2  r~ 5 - 2 
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exp 
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(22) 
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(Note  that  the  units  of  I(t)  are  energy  (pulses)  per  unit  time  and  of  G(t)  , 
energy  squared  per  unit  time.) 

Suppose  the  maximum  and  minimum  channel  contents  near  the  origin  (t  *  0) 
are  G^^  and  G  -  these  correspond  to  values  of  the  cosine  term  of  approxi¬ 
mately  +1  and  -1  respectively.  Then  if  t  is  the  sample  time,  ie  the 
width  of  a  channel,  and  if  m  is  taken  to  be  unity. 


^  2  /ir  r0 

C*  2  /u2  +  v2 

from  which  c^  can  be  found.  Also 

G(0) 

giving  c2  . 

The  total 
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2  3/tt 
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2  4  2 
V  u  +  V 
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number  of  pulses  received  by  the  correlator  is 
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V2  rr~ 


-T/2 


1  V  2  rr—2  +  c2T 

/  u^  +  v^ 


Hence  the  number  of  signal  counts  is 


N 

s 


r0 


/72  +  v^ 


(23) 


and  the  number  of  background  counts  is 


Nb  =  C2T 


In  Ref  10  a  correlation  function  acquired  over  a  period  of  25  ms  is  dis¬ 
played.  For  this  experiment,  the  relevant  values  are  found  to  be,  approximately, 

Gmax  *  22»  Gmin  *  1 1  •  The  sample  time  setting  t  was  0.25  us  ,  the  beam 
•  ®  .1 

radius  about  500  microns  and  the  velocity  across  the  beam  about  33  m  s  .  Wte 

find  for  this  case 

c.  =  l .819  us  ' 

-1 

c2  =  0.0283  us 

N  =  34 

s 

Nb  -  708 
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The  peak  count  rate  would  have  been,  putting  t  =  0  and  m  =  1  in  (21) ,  3.67  MHz, 
while  the  background  count  rate  was  28.3  kHz.  It  is  interesting  to  note  that  a 
good  velocity  estimate  appears  to  be  obtainable  from  only  34  signal  pulses,  in 
the  presence  of  considerable  noise. 

4  CONCLUDING  REMARKS 

It  has  been  demonstrated  that  the  characteristics  of  the  measuring  volume 
make  some  degree  of  approximation  in  the  interpretation  of  the  Doppler-difference 
signal  or  of  its  autocorrelation  function  unavoidable,  and  that,  in  the  absence 
of  a  complete  model  for  the  velocity  field  of  the  fluid  under  study,  further 
simplifying  assumptions  have  to  be  made.  The  errors  which  result  can  be  made 
small,  however,  by  careful  design  of  the  experiment.  Specifically,  the  optical 
geometry  should  be  arranged  wherever  possible  so  that  the  period  of  a  Doppler- 
difference  cycle  is  very  much  less  than  the  beam  radius  transit  time;  this  case 
involves  the  smallest  number  of  assumptions  about  the  nature  of  the  flow.  In 
addition,  it  has  been  shown  that  the  data-reduction  procedure  is  then  of  the 
simplest  kind. 
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